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Abstract: We propose a direct and real-time ranging scheme using 
an optical frequency combs, able to compensate optically for index of 
refraction variations due to atmospheric parameters. This scheme could be 
useful for applications requiring stringent precision over a long distance 
in air, a situation where dispersion becomes the main limitation. The 
key ingredient is the use of a mode-locked laser as a precise source for 
multi-wavelength interferometry in a homodyne detection scheme. By 
shaping temporally the local oscillator, one can directly access the desired 
parameter (distance) while being insensitive to fluctuations induced by 
parameters of the environment such as pressure, temperature, humidity and 
CO2 content. 
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The accuracy of precise length measurements is commonly limited by dispersive effects. For 
instance, the dispersion of air is a crucial issue for geodetic surveying |T| or for the optical 
link between a ground station and a satellite [2J. The lack of knowledge about atmospheric 
parameters can then act as the main limitation to optical length measurements. 

Numerous groups around the world tackle this issue of long range distance measurement. 
For instance, in atmospheric links such as satellite ranging or Lunar Laser Range, the accuracy 
provided by time-of-flight measurements reaches the millimeter level or below J2][3). Interfer- 
ometric distance measurements lead to potentially increased accuracies. However, for simple 
interferometric experiments, the periodicity of the signal results is ambiguous up to an absolute 
distance of a half-wavelength. This ambiguity distance can be extended by combining signals 
from multiple wavelengths. For instance, the Very Large Telescope Interferometer uses dual- 
field interferometry to reach nanometric accuracy in a 100 m air- filled delay line (0) . 

It is also possible to combine both time-of-flight and phase measurement to obtain a better 
sensitivity and an absolute distance measurement. The ideal tool for this is the phase-stabilized 
mode-locked femtosecond laser which delivers a frequency comb that can be seen as a perfect 
source for multi-wavelength interferometry |5]|6]|2l. The problem we address here is the way 
to use this tool in a complex environment, in order to perform a dispersion free measurement. 
We will more precisely treat the example of distance measurement in air independent of the 
variation of physical parameters such as pressure, temperature, CO2 content or humidity. 

In a multi-parameter environment, where many physical factors can affect the accuracy of 
distance measurements, these extra parameters need to be measured and their effects com- 



pensated. This is the general strategy of the well-known multicolor schemes |8|, which we 
introduce in the first part of this paper. 

In the second part, we derive fundamental limits for optimal measurement schemes (i.e. 
reaching the Cramer-Rao bound with coherent states (9J \W[ \TTl ) for distance measurements 
through a dispersive medium, using mode-locked lasers. The technique is based on temporal 
mode-dependent interferometry. We show that, in contrast to multicolor schemes for instance, 
only one measurement is necessary whatever the number of parasitic parameters we want to 
cancel. This is done at the cost of a precise spectral mode shaping of the frequency comb that 
is used. 

We finally propose a general all-optical experimental setup using pulse shaping and homo- 
dyne detections to reach the previous limit. Our distance measurement can be made insensitive 
of the environment parameters such as temperature or pressure. No post-processing is needed to 
achieve a measurement limited by the laser noise. Another advantage over existing schemes is 
that mode-locked lasers are remarkable tools for optical measurement because of their intrinsic 
high stability. 

Distance measurement in air : multicolor schemes 

Let us start by describing the general technique of multicolor schemes. If one needs to mea- 
sure a distance in air, the fluctuations of parameters such as pressure or temperature limit the 
achievable accuracy. The general idea to solve this problem is to perform several measurements 
at different wavelengths to gain informations about these parameters and compensate the meas- 
ured length for their variations. For example, in the two-wavelength interferometry (2WI), a 
given distance L is measured using two wavelengths X\ and A 2 - The two observables and 
deduced from the measurement are such that ll8l [T2l[T3ll 

Lfa = n${X\)L and L fc =«0(A 2 )i- (1) 
For dry air (pressure of water vapor P w — 0), one finds 

L = L^+a(L h -L$ 2 ) with a = J K ^i) ' (2) 

where K(X) is given in Appendix|A] and is calculated from the Edlen model of air. The param- 
eter a is independent from pressure, temperature and CO2 content of air. The standard quantum 
limit for a distance measurement is obtained by taking (dL^) — 2 J^ m ■ The sensitivity 8L 
on the distance then strongly depends on the value of a. Typically, for Ai = 1064 jim and 
A2 = 532 fim, a is of the order of 60. Hence, for a mean photon number iV = 4x 10 16 per 
wavelength (10 mW power and an integration time of 1 s), one gets 

(5L)| h ^~3x 10- 14 m. (3) 

This is the shot noise limit in distance variation measurement (displacement), but one should 
note that for an absolute distance characterization (ranging), the precise knowledge of a is the 
main limitation factor (8a/ a ~ 1%). 

For moist air, one can no longer use Eq. ©. If the pressure of water vapor is unknown or 
changing over time, it leads to a systematic error (see lfT3l for a discussion). A solution is then 
to consider a third wavelength A3 lfT4l and a third measurement so that 

L = L <h +P{L <h -L^) + Y(L h -L <h ). (4) 

Expressions of /3 and 7 can be found in lfl4l and do not depend on pressure, temperature, 
CO2 content and humidity. Here again these factors can be large, reducing the sensitivity of 



displacement measurement. For a same total number of photons and typical wavelengths, the 
shot noise limit is now around 10~ 12 m. In addition to this degradation, the three-wavelength 
scheme is experimentally more involved. 

From this example one sees that both the required sensitivity and physical characteristics of 
the medium conditions on the number of extra parameters one has to take into account. For 
each of these parameters an additional measurement at a different wavelength is necessary. 

1. Efficient measurement through dispersive media 

We will now give a more general and systematic approach to this problem. First, we derive gen- 
eral equations for the propagation of light through a dispersive medium whose characteristics 
depend on external parameters. We then elaborate on the ideas developed in J6] [TT] and deal 
with a very general approach on how to efficiently measure parameters affecting the propaga- 
tion of a light pulse. We derive fondamental limits imposed by the quantum nature of light. One 
should note nevertheless that we limit ourselves to the study of coherent states, non-classical 
quantum states being beyond the scope of this article. 

1.1. Propagation in a dispersive medium 

We consider the propagation of an electromagnetic field along the z direction in a weakly dis- 
persive medium. Its propagation from a source to a detector can be affected by a given set of 
parameters p — {p i , . . . , pt , ... ) that modify the propagation distance L and/or the characteristics 
of the dispersive medium: these may be environmental parameters such as air pressure, tem- 
perature, etc., or a physical displacement of the source (or the detector). Neglecting polarization 
effects and using the paraxial approximation, we assume the field to be in a single transverse 
mode (such as a TEMoo mode), and thus will not write the transverse dependence of the field. 
We further assume Fourier-limited pulses (assuming perfect temporal coherence) and write the 
dispersed field, as seen by the detector, as a scalar field: 

*(t,p) = 4>u(t,p), (5) 

where u(t,p) is the normalized mean field mode (integrated over the measurement time of the 
detector) and Sq is a normalization constant which depends on the number of photons N. 
In the following, it will be more convenient to work in the Fourier space: 

g(co,p) = j g(t 1 p)e im dt, u(co,p) = j u{t,p)e ia "dt . (6) 

We define C0q and Acq as the mean value and variance of the field: 

do = J co\u(co,p)\ 2 dco, Aco 2 = J "(to — coo) 2 \u(co,p)\ 2 dco . (7) 

Let us consider an input field $i(co) whose frequency profile is known. For the sake of sim- 
plicity, the field is considered Gaussian (the same final result can be reached with a non Gaus- 
sian field but it involves more complex calculations). This field propagates on a distance L(p) 
through a dispersive medium with a dispersion relation k(co,p) which depends on p through the 
refractive index n§ ((0,p). In the frequency space, this propagation is characterized by a spectral 
phase k(co,p)L(p): 

S(co, P ) = ^o))e'^™ k(co,p) = n ^ (0 ^ (0 . (8 ) 

c 

In this paper, we neglect any absorption of the medium, i.e. we consider a real refractive index. 
The previous equation ([H} is therefore also valid by replacing fields § by normalized modes u. 



For a weakly dispersive medium, the dispersion relation k(co,p) can be expanded to the 
second order around the mean frequency COq: 



(Oq (p) + (CO - fflo) f g (p) + ^ J 3 ^ f G VD (P) 
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n g is the group index, and 2 \n'^{(OQ,p) + ^/^(filo, j?)J corresponds to the Group Velocity 
Dispersion (GVD). 

In the temporal domain, the field £{t,p) — e~ l(CK>t §(t ,p), the envelope § of the field travels 
at the group velocity ~ while the carrier o>o moves at the phase velocity ; a non zero group 
velocity dispersion leads to a broadening of the envelope. 

Any change of a parameter in p that affects the distance L(p) will contaminate all quantities 
t$, t g and ?gvd, as well as any variation of the refractive index of the medium. In Section[2l we 
show how to uncouple, in air, variation of L from variation due to four different environmental 
parameters: pressure, temperature, humidity and CO2 content. Note that a generalization to 
other environmental parameters can be obtained by expanding equation (O to higher orders of 
the spectral phase and applying the methods developped later in the paper. Nevertheless, these 
4 parameters are typically the most relevant for air. 

1.2. Detection scheme and Cramer-Rao bound 

The general problem of estimating a parameter pi E p encoded in a light beam S"(p) has been 
treated in J9][l5j. The ultimate limit of sensitivity in the measurement of p\ is given by the so- 
called quantum Cramer-Rao bound, which can be computed once we specify the quantum state 
of the light beam (coherent state, squeezed state, entangled state etc..) ifTTI . For Gaussian states, 
this Cramer-Rao bound can be reached experimentally with a balanced homodyne detection 
scheme, as represented in Figure Q] The general idea is that the homodyne detection signal is 
proportional to the projection of the input field into the Local Oscillator (LO) mode. 
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Fig. 1. Detection scheme for measuring p, at the Cramer-Rao bound. 



For a small variation of the set of parameters p, the field reads: 



«#(o),0^/^<f(^ , (13) 



where Wj are normalized modes such as w,-((0) = jrj^((O,p = 0), and K~i are dimensional 
normalization constants. Introducing the standard L 2 inner product (f,g) = J f* (co)g(co)da), 
one simply has K, = yj(^, g). 

One should note that in general the modes w, do not form an orthogonal basis. This implies 
that independent measurements of each parameter become a complex problem that we now 
discuss in detail extensively. Let us first consider the case where only one parameter pi € p is 
influencing the length measurement. It is shown in lfl6ll that in a homodyne detection scheme, 
if the LO mode is proportional to w>j (and if there is no phase difference between the LO and 
the signal <S), the detected signal S[w,] is given by 

S[ Wi ] = ^Re[(u(p), Wi )} = Pi . (14) 

For a coherent state illumination with a mean photon number N, the noise in the measurement 
is Api = 2y/NKj. Therefore, the smallest pi that can be measured (i.e. for a signal to noise ratio 
equal to one, S/Ap = 1) is given by 

(/,i)min= ^W (15) 

This value coincides with the Cramer-Rao bound with coherent states IfTTII . This shows that a 
homodyne detection with a LO shaped in mode w, defines what is called an efficient measure- 
ment of pi. Moreover, one sees from the previous expression that the sensitivity of the measure- 
ment depends both on the number of photons N and on K{, the latter reflecting the characteristic 
variation of the mode with the parameter p,-. From now on, the mode wi will be called the 
detection mode of the parameter p,-. 

An experimental demonstration of the efficiency of such a scheme has been realized for pa- 
rameters corresponding to spatial displacement of a beam fl7l[l6l . and a theoretical proposition 
for a time delay through a dispersion-less medium has been made in 0. One should stress that, 
generally speaking, these kind of experiments are sensitive to variation of parameters within the 
detection bandwidth, limited here by light time travel. On the other, this system is immune to 
any fluctuations, whatever their frequencies, of parameter corresponding to orthogonal modes, 
given we stay in the linear regime. 

Let us now consider the general case where there exists at least another parameter p\f\ such 
as Wi and Wj are not orthogonal. In that case, a homodyne detection with LO mode w, will also 
be sensitive to pj. We show here that this issue can be resolved by 'purifying' the detection 
mode Wi into a new mode wf that is now orthogonal to Wj. This new shape of the LO allows to 
measure p\ independently of pj. Nevertheless, because it differs from the detection mode w,, 
it leads to a reduced sensitivity in the measurement of pi. In the general situation, the purified 
mode for a given parameter p; is orthogonal to the hyperplane formed by {w,%£;}, and the 
normalization factors are given by Kj 1 = Ki(m,Wi) < Kj. The sensitivity in the measurement of 
Pi is therefore decreased and given by (p;)mi n = 2 >/Nk!' ' 

The choice for the mode of the LO is an experimental tradeoff between accuracy (no per- 
turbation coming from pj) and precision. We further elaborate on this point in the following, 
taking as a simple example the measurement of u, % and ?gvd introduced previously. 

To this aim, we introduce a controlled perturbation p^ <C of the phase delay f<j, — > t^+p^, 
a perturbation p g -C t g of the group delay t g — > t g + p g and a perturbation pgvd "C ?gvd of the 
group velocity dispersion delay ?gvd ? gvd + Pgvd- The corresponding detection modes are 



given by: 



w^(fi)) = iu(a>) = vo(ffl) (16) 

w g (fi>) = i Am M(fl))=vi(fl)) (17) 
/ x 1 ((o-a>o) 2 / \ 1 . , fl , . 

WGVD(fi)) = l-j= Ao}2 H(ffl) = -^=Vo((H) + y-V2(fi)) (18) 



where we have introduced {v,(a))} the orthonormal basis of spectral Hermite-Gaussian modes 
whose expressions are given in Appendix IB] The normalization factors are given by Ka = (Oq, 

K g =A(0 and/T G VD-V3^. 

It is clear that w$ and wgvd are not orthogonal, which implies that a measurement using a 
LO Wtj, will not be accurate because of its sensitivity to pgvd, and vice versa. More precisely, 
measurements over the various detection modes will yield the following signals: 

S[wq] = P<f, + -^PGVD (19) 



S[w g ] = p g (20) 

1 CO 2 

S[w G \i D ] = 2~s^2P<l>+PGVD- (21) 



In order to measure only p$, one has to consider the purified mode introduced previsously, 
which is orthogonal to the detection modes of the other parameters. Formally, w^(co) is or- 
thogonal to the hyperplane generated by {w g , wgvd} in the vector space {w^ , w g , wgvd} and is 
given in the present case (up to a scalar factor) by: 

^(») = y|v (£o)--^v 2 (a)). (22) 
A measurement with LO wj yields: 

S[w$]=Pt. (23) 

2 



Defining — (w^,wq) — y |©o, the sensitivity of this measurement is: 



The two previous equations show that it is possible to retrieve a phase delay information 
insensitive to any group velocity dispersion fluctuations with only one homodyne measurement. 
This improvement in accuracy is made at the cost of a decreased precision, determined by the 
overlap between the purified and non purified modes (here the degradation is given by W3/2). 

The same process can be applied to obtain the purified mode for measuring pgvd- 

<v D ( fi) ) = v '2(«) with ^gvd = ^2— ■ (25) 

(Oq 

The link between the detection modes and the purified modes is shown in Figure |2] and in 
TableQ] 




2. Application to the measurement of a distance in air 

Let us now consider the specific case of measuring a distance in air, independently of the fluc- 
tuation of environmental parameters such as pressure and temperature. Indeed, induced index 
refraction fluctuations are the main limitations to precise distance measurement. To access the 
absolute length, one needs to know precisely the air index variation with these parameters. One 
solution is to measure precisely these parameters and use an air model to calculate the refractive 
index, for instance the Edlen [18] or the Ciddor [19| equations. But these techniques are not 
immune to local variation of the parameters, and some parameters such as the partial pressure 
of water vapor are very difficult to precisely access. One can also measure directly the local 
refractive index using a refractometer ll20l . 

Another solution is to make several measurements in order to compensate for these effects. 
This is the principle of multicolor interferometry |8 |, whose state of the art is based on second 
harmonic generation ll2Tl[T2ll . as we have presented in the introduction. Here we develop our 
new experimental scheme based on mode-locked lasers interferometry as introduced in the 
previous section. This technique is another kind of multicolor interferometry and allows for 
direct measurement of distance in air independently of parameters from the environment. 
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Table 1. Summary of the different LO modes and the sensitivities. 



2.1. Ranging through air 

We apply the technique developed in Section [L2l to the measurement of the absolute distance 
L in air. The parameter to be measured with high sensitivity is pi = L. Fluctuations of the 
environment do perturb this measurement. They can be separated into two groups of parameters. 
Firstly temperature T, pressure P and CO2 content x affect air index though the same function, 
as can be seen in the air model developed in Appendix lAl They will be described by only one 
parameter px = X(T,P,x). Secondly, pressure of water vapor P w has an independent influence, 
for which we define the parameter pp w = P w , 

One can calculate the corresponding detection modes using the Edlen model of refractive 
index recalled in the appendix and the second order development of the electric field introduced 
in the first section. The distance detection mode is given by : 

w L {co) = — (<Hovo(o)) + Aa)vi(£o)) (26) 

and the two other detection modes are given in the appendix. 

These modes are not linearly independent. Thus if we consider a possible experimental 
scheme with homodyne detection in the detection mode for L (see Figure O the measured 
signal will be : 

Kx Kp 
S[w L ] = pl + —{wl,w x )px + -tt-(wl,wp w )pp w - (27) 

Therefore, the signal will be contaminated by variations of the different parameters px and pp w . 
Let us first compute the shot noise limit in the case where px and pp w are zero (or sufficiently 
small). To compare with multicolor schemes, in the remainder of this section the measure- 
ment is performed using = 8 x 10 16 photons and assuming a laser bandwidth of Aco = ^p- 
(corresponding to 3 fs FWHM Fourier-limited pulses). The shot noise limited sensitivity to 
displacement is about 2 x 10~ 16 m, comparable to usual interferometric measurement schemes. 
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Fig. 3. Direct distance measurement with an appropriately shaped LO. 



One can evaluate the contamination from the other parameters calculating the pre-factor of 
px and pp w in equation d27l ). One finds ^j^{wl,Wx) — 27x 10~ 5 Pa -1 and j;j^-{wl,Wp w ) = 
—3.7 x 10 _10 Pa _1 . These factors are big compared to shot noise limited pure distance measure- 
ments. It means it is necessary in this case to take into account the parasitic parameters in order 
to preserve the accuracy of the measurement. To solve this issue, we can apply the detection 
mode purification procedure introduced previously : 

P(„\ ir,\ {wl,wp w )-{wx,wp w ){wl,wx) , , 

oewi(fi)) - -, To W Pw {CO) (28) 



{W X ,Wp w 

( WL ) wx ) " (wx , wp w )(wl, w Pw ) 
1 - (w Xl w Pn ,) 2 



wx((o) . (29) 



The spectral profiles of the purified and non-purified modes are plotted in Figure [4] The nor- 
malization constant as well as the derivation of this mode can be found in Appendix [C] It is 
found that : 



k p =Kl i ( w l,w Pk ) 2 + (wl, wx} 2 - 2(w L , w Pw }(w L ,wx)(wx , w P J 
L L \j l-(w x ,w Pw ) 2 

In that case, the shot noise limit has the following value : 

(5L)f H °^- 7 |-=2xlO- 11 m. (31) 

As a matter of comparison, a purification for px only (which would be equivalent to the two- 
color scheme) leads to a shot noise limit of (5L)j^ = 2 s/J/k p = ^ x ^ 13m - Even if the 
precision is 2 orders of magnitude better, the accuracy of the measurement is not controlled 
because of the unknown value of humidity introducing a systematic error. 

Using that scheme, one can perform a single-shot measurement of distance immune from 
environmental parameters fluctuations. Sensitivity is of the same order of magnitude than in 
multicolor scheme (slightly lower in our examples, but this is simply due to smaller spectral 
width), and depends on how many parameters one wants to get immune from. From a more 



Fig. 4. Spectral profile of modes wi and w p L for 3 fs FWHM Gaussian pulses. 



technical point of view, no complex frequency generation is necessary, a femtosecond source 
is enough to perform the measurement. Furthermore, the scheme described here could be ex- 
tended to more parameters. Of course, the mode shaping can be more difficult to produce and 
in particular the precision required in shaping becomes more stringent when more parameters 
are considered. In a realistic implementation, one can set up active pulse-shaping with search 
algorithms to determine the proper shaping. It is significant to notice that all the work is to be 
done at the detection stage and not on the light sent through the medium, which makes it much 
easier to handle. 

3. Conclusion 

We have exhibited a novel optimal and all-optical scheme to measure in real time a distance 
compensated for refractive index changes. It relies on a homodyne detection whose local oscil- 
lator projects the measurement on an appropriate mode, hence no post-processing is necessary. 
We believe this is a simplification compared to existing schemes such as spectral interferometry 
for example, where derivatives of the spectral phase have to be done after the measurement. 

Let us finally mention that this scheme can further be improved in order to go below the 
standard quantum limit. It is well known that the sensitivity of a measurement can outreach 
the shot noise limit by using quantum resources such as squeezed or entangled light ll22l . In 
the scheme presented in this paper, this can be achieved by using a multimode signal light 
beam with squeezing in the detection mode associated to the measurement, as demonstrated in 
Ref. JED. 

A. Index of refraction of air 

The common equations used to derive the wavelength dependence of the refractive index of air 
are given by the Edlen equations lfT8l . modified by different authors since that time [19 23ll24l 
|25] |26] |20) . The accuracy of those equations are roughly of the order of a few 10~ 9 for dry air 
and 10~ 8 for moist air and needs to take into account a large number of parameters, usually 
temperature, pressure, CO2 content, pressure of water vapor. Moreover, different studies do not 
necessarily agree and do not cover the whole spectrum. A recent precise measurement of the 
refractive index of air has been done in [27 1 . Experimentally, measuring all those parameters 
can be difficult for certain situations, and in addition those parameters have to be measured in 
real time, in order to compensate for fluctuations of the refractive index. 



In this paper we consider the updated Edlen formula of Bonsch and Potulski |20| : 

n$(o,T,P,x,P w ) - 1 = K{a)X(T,P,x)-g{a)P w (32) 
where a = 1 / A is the wavenumber and 
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K ^ = 10 l A+ 130^ + 38^J ' (33) 

P I + 10~ & (E — F T)P 
X(T,P,x) = j-^- , -[l+H(x-0.M%)], (34) 

g(a) = 10- 10 (/-/cj 2 ) , (35) 

with a in jtim -1 , P and P w in Pascal (Pa), T in degree Celsius (°C) and x, the CO2 content, in 
percentage. The different coefficients further read 

A = 8091.37 , B = 2333983 , C= 15518, (36) 

D = 932164.60 , £ = 0.5953 , F = 0.009876 , G = 0.0036610, (37) 

H = 0.5327, (38) 

7 = 3.802 , 7 = 0.0384. (39) 

The group index is therefore given by 

n g (a,T,P,x,P w )-l = {K(a) + aK'(a))X(T,P,x)-(g(a) + ag'(a))P w , (40) 

so that n g -n$ = o {K'{a)X(T,P,x) - g'(a)P w ). 

B. Hermite-Gaussian set of spectral modes 

For a Gaussian mean field mode 

1 1 ("-"o) 2 



/Aco (2 s; 

we introduce the Hermite-Gauss modes 

v n ((0) = i^L=H n ) u(co). (42) 



V \/2Aco 

These modes {v„(co)} form an orthonormal basis of modes. 

In order to describe the spectral phase up to the k th order, one needs to use the basis {v„(co)} 
up to the mode k. In this paper, we develop the spectral phase to the second order; therefore, we 
use the following modes: 

vo(ffl) = iu((0) (43) 

vi (ffl) = i^-^u(co) (44) 
Acq 

»<•> - ^(t-'W 

The mode spectral profiles are plotted in Figure [5] 



Fig. 5. Spectral profile of Hermite-Gaussian modes vq, vi and V2. 



C. Detection modes of environmental parameters 

Detection modes read : 
1 

~ck L 

LK(coq) 



W L (O0) = 

w x {co) = 

Wp w (<X>) = 



(coovo(co) + Afi)vi(fi))) 



(46) 



cK x 
cK Pw 



fflo + — (Si + 82)) v (o))+A(o(l + 5i)vi((o) + \/2— (Si + 52)v 2 ((o) 
oh J coq 

£O0 + ^-(TJl + TJ 2 ) N ) V (£»)+A£0(l+T}i)vi(a)) + V2^-(TJi+T} 2 )V2((») 



where we have defined characteristic quantities that do not depend on the environmental pa- 
rameters T, P, x and P w : 
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Defining : 

K L = 

K x = 
Kp„ = 



1 
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K(ooq)L 
c 

8(ooo)L 



O0q + Aoo 2 



(47) 
(48) 
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Aoo 2 
fflo 



2 a 4 

(5! + 52)) +Aoo 2 (l + 8 l ) 2 + 2-^ r (8 l + 8 2 ) 2 (50) 

fflT 
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Aoo 2 
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Aoo 4 . 



(77! + 77 2 ) +A0) 2 (1 + T7 1 ) 2 + 2^(77! + T72) 2 (51) 



cat 



Measurements with the detection modes give: 

K X 



M[wl] 
M[w x ] 
M[wp w ] 



Kp w 



PL + —{WL,W X )px + —{WL, W Pw )pp w 

— (vf£ , w x )pl + Px + -7- (wl , wp w ) pp w 

Kx Ax 



K x 



{WL,Wp w )p L -+-—{WX,Wp n .)p X +pp lr 
Kp w Kp w 



(52) 
(53) 
(54) 



Calculation of the purified mode w p L : we first calculate a orthonormal basis in the plane 
{wx , wp K }, which leads for example to a basis {wx , w' Pw } where 

w'pJco) = = (w Pw ((0) - (w x ,wpjw x ) (55) 

^l-(w x ,wpj 2 

Then we do a Gram-Schmidt orthonormalization of {wx , w' P , wl} which gives: 

wl(co) cx w L (co) - (w L ,wx)w x ((0) - (w L ,w l pJwp w (56) 
The normalization constant is: 



( wl , w Pw ) 2 + ( w L , w x ) 2 - 2 (w L , w Pw ) (w L , w x ) ( w x , w Pw ) 



1 - (w Xi wp w ) 2 
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